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ABSTRACT 

We present numerical simulations designed to test some of the hypotheses and 
predictions of recent models of star formation. We consider a set of three numerical 
simulations of randomly driven, isothermal, non-magnetic, self-gravitating turbulence 
with different rms Mach numbers Ms and physical sizes L, but with approximately the 
same value of the virial parameter, a « 1.2. We obtain the following results: a) We 
test the hypothesis that the collapsing centers originate from locally Jeans-unstable 
("super- Jeans"), subsonic fragments; we find no such structures in our simulations, 
suggesting that collapsing centers can arise also from regions that have supersonic 
velocity dispersions but are nevertheless gravitationally unstable, b) We find that 
the fraction of small-scale super- Jeans structures is larger in the presence of self- 
gravity, c) There exists a trend towards more negative values of the velocity field's 
mean divergence in regions with higher densities, implying the presence of organized 
inflow motions within the structures analysed, d) The density probability density 
function (PDF) deviates from a lognormal in the presence of self-gravity, developing 
an approximate power-law high-density tail, in agreement with previous results, e) 
Turbulence alone in the large-scale simulation (L = 9 pc) does not produce regions 
with the same size and mean density as those of the small-scale simulation {L — 1 
pc). Items (b)-(e) suggest that self-gravity is not only involved in causing the collapse 
of Jeans-unstable density fluctuations produced by the turbulence, but also in their 
formation. 

We then measure the "star formation rate per free-fall time" , SFRff , as a function 
of Ms for the three runs, and compare with the predictions of recent semi-analytical 
models. We find marginal agreement to within the uncertainties of the measurements. 
However, within the L = 9 pc simulation, subrcgions with similar density and size to 
those of the L = 1 pc simulation differ qualitatively from the latter in that they exhibit 
a global convergence of the velocity field V • t> ~ —0.6 km s~^ pc~^ on average. This 
suggests that the assumption that turbulence in clouds and clumps is purely random 
is incomplete. We conclude that a) part of the observed velocity dispersion in clumps 
must arise from clump-scale inwards motions, even in driven-turbulence situations, 
and b) analytical models of clump and star formation need to take into account this 
dynamical connection with the external flow and the fact that, in the presence of 
self-gravity, the density PDF may deviate from a lognormal. 

Key words: interstellar matter - stars: formation - turbulence 



* E-mail: e.vazquez@astrosmo.unam.mx 
t E-mail: rf.gonzalezOastrosmo.unam.mx 
X E-mail: j.ballesteros@astrosmo.unam.mx 
§ E-mail: a.gazol@astrosmo.unam.mx 
^ E-mail: jskim@kasi.re.kr 



1 INTRODUCTION 

The role of the velocity field in the process of star 
formation (SF) remains not fully understood. It is 
generally believed nowadays that supersonic turbulent 
motions of a given size scale in molecular clouds (MCs) 
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shoul d ha ve a dual role in relation to SF (e.g.. ISasao 
19731 : iBallesteros-Paredes. Vazguez-Semadeni. fc Scale 



19991: I Vaz guez-Semadeni fc PassotI 

Klessen. Heitsch fc Mac Low 2000i : iMcKee fc Ostrikeij 



19991: 



20071 ): they are thought to provide support towards scales 
larger than their own, while simultaneously promoting 
collapse of smaller scales. However, the origin and nature 
of the observed supersonic mot ions remains contro v ersial . 
More than three decades ago, iGoldreich fc KwanI (119741 ) 
suggested that the observed supersonic widths of various 
molecular lines could be representative of gravitational 
contraction, al though this suggestion was di smissed shortly 
thereafter by IZuckerman fc Palmen (|l974l ) through the 
argument that if all MCs were collapsing and converting 
most of their mass into stars in roughly one free-fall time, 
the resulting star formation rate would be at least 10 
times larger than that pre sently observed in the Galaxy. 
IZuckerman fc EvansI l| 19741 ) then proposed that the super- 
sonic linewidths in the clouds are produced primarily by 
local, small-scale motions, a scenario to which we refer as 
"local turbulence", and which has been widely accepted 
until recently. Built into this scenario is the notion that 
local turbulence acts as a sort of isotropic pressure, so that 
it provides an important (and perhaps dominant) contribu- 
tion in the support of the clouds against their self-gravity, 
in particular allo wing for a state o f near hydrosta tic equi- 
librium (e.g.. McKee fc TanI l20q | : | Hennebelle fc Chabrie3 
l2008l : lField. Blackman. fc Ketoll2008l V 

However, this scenario f aces a numb er of problems. 
First, it is well known (e.g., I FrischI 119951 ) that turbulent 
flows possess the largest velocity differences at the largest 
scales, a property which implies that the largest velocity 
differences within a turbulent structure are expected to 
occur at scales comparable to the size of the structure 
itself. That is, Lagrangian (i.e., moving with the flow) 
clumps are expected to be continually distorted (sheared 
and/or compressed) by these large-scale motions, in 
contradict ion with the idea of their being in near equi- 
libriu m (JBallesteros-Paredes. Vazguez-Semadeni. fc Scald 
I1999I ). Such large-scale motions are indeed observed in 
molecular clouds and their s ubstructure (e.g.. iBrung |2002| : 
lOssenkopf fc Mac Lowll2002l : iHever fc Bruntll2007h . 

Second, there is the problem of how to main- 
tain the observed turbulence levels. Early suggestions 
were that the tur bulent motions consisted of hydromag- 
netic waves (e.g.. lArons fc MaxI Il975l : iMouschoviasI 1 19761 : 
IShu, Adams fc Lizand [l987|), which would be less dissi- 
pative than supersonic hydrodynamic turbulence. How- 
ever, numerical experiments demonstrated that MHD tur- 
bulence generally decays as fast as hy d rodynamic tur- 
bulence Istone. Ostriker fc Gammi3 Il998l ; iMac Low! 1 19991 : 
iPadoan fc Nordlundl Il999l ). except perhaps if it is un- 
balanced (i.e., if the energy flux along field lines in 
one direction differs from that in the opposite direction; 
ICho. Lazarian fc Vishniad 120021 ) , although it is not clear 
whether this applies to MCs. 

Maintenance of the turbulence in MCs by stellar energy 
feedback has also been proposed, although it is not yet clear 
whet her this feedback can keep the clouds near equilibrium 



2OOOI : 

2007! ) 



Norman fc Silklll980l:[McKee"l989':lMatzner fc McKed 



Krumholz. Matzner. fc Mc Kcc 2006; Nakam ura fc Lil 
20071 ) or rather disrupt them (e.g., IWhitworthI 1 19791 : 



LarsonI [l987l: iFranco. Shore fc Tenorio-Tarid Il994l : 
Hartmann. Ballesteros-Paredes fc BerginI I2001I ). In both 



cases, the stellar feedback has been proposed as a regu- 
lating mechanism of the star formation efficiency (SEE). 
In the former, because the high stationary turbulence 
levels are expected to maint ain the star formation rate 
(SFR ) at low values (e.g., Klessen. Heitsch fc Mac Lowl 
2OOOI: iHeitsch. Mac Low fc KlessenI I2OOII 



Vazguez-Semadeni et al.l |2003| . l2005bl : iNakamura fc Lil 



20071 ). while in the latter the SFR can be large after a cloud 



forms, but the cloud is soon dispersed after SF begins. In- 
deed, clusters over 5-10 Myr o ld are generally observed to be 



devoi d o f molecular gas (e.g., Leisawitz. Bash, fc Thaddeua 



19891: Hartmann. Ballesteros-Paredes fc BerginI I2OOI 



Ballesteros-Paredes fc Hartmannll2007^ . 



Another alternative that has been proposed recently is 
that MCs obtain their turbulence from the compressions 
that fo rm them in the first place through various insta- 
bilitie s iJKovama fc In utsuk all2002 ;_ Vazguez-Semadeni et al.l 



Audit fc HenncbcUe 2 0051 : Heitsch et all |2005| , I2OO6I : 



IVazguez-Semadeni et al.l I200d ) , although in this case the 
turbulence should begin to decay after the compression that 
formed the cloud subsides. IVazguez-Semadeni et al.l (|2007| ) 
showed simulations in which, when this happens, the cloud 
begins to contract gravitationally, exhibiting a virial-like en- 



ergy balance | Eg 



2i5kin while doing so, even though the 



cloud is never in virial equilibrium. This can in fact be un- 
derstood because, in the case of gravitational collapse, the 
velocity dispersion is at most within a factor of V2 from 
the value needed for equipartition. Thus, both a system in 
virial equilibrium and one in free-fall are, to order of magni- 
tude, in energy equipartition, and are, thus, observationally 
indistinguishable on the basis of velocity dispersion alone 
dBallesteros-Pareded IJOOJI ) . 

If the simulations by IVazquez-Semadeni et al.l (|2007l ) 
are representative of actual MCs, the nonthermal motions 
in the clouds could transit from being initially due to actual 
random turbulence, to later being due to gravitational con- 
traction. In this case, the pseudo- virial energy balance would 
be a manifestation of the gravitational contraction rather 
than of virial equilibrium. We refer to this mode of fiow, in 
which the velocity field at all scales contains a significant in- 
flow component, as large-scale inflow (LSI). Such a flow can 
be due to either generic dynamic compressions in the ISM 
(e.g., expanding HII regions or supernova remnants, or the 
transonic turbulence in the warm ISM), or to various large- 
scale i nstabilities, such as gravitati onal or magneto- Jeans 
(e.g., JLi. Mac Low fc Kl essen 20051: iKiin fc Qstrikeill200d . 
see also the review by Hennebelle, Mac Low fc Vazquez- 
Semadeni 2008). 

The scenario of LSI driven by self-gravity is actually 
frequently encountered. Since it is standard in simulations 
of star formation in clouds with decaying turbulence (e.g., 
JBate. Bonnell. fc Brommll2003l : JBonnell fc Batdl200d ). it is 
often associated with a regime of decaying turbulence, al- 
though in principle there is no reason why it should be 
only applicable in this case. JField. Blackma n. fc Ketd (120081 ) 
have attempted to describe a gravitationally-driven mass 
cascade that involves the formation of smaller-scale struc- 
tures by gravitational contraction of larger-scale ones, fol- 
lowed by virialization at the small scales. The LSI scenario 
is also consistent with a number of recent studies suggest- 
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ing that indeed some MCs IIHartmann fc Burkertll2006l ) and 
clumps (|Peretto. Hennebelle fc AndrgboOTI ) may be under- 
going global gravitational contraction. It should be noted, 
however, that the simulations supporting this scenario have 
generally been n on-magnetic, possibly biasing the results. 
lElmegreenI (|2007l ) has recently suggested that clouds may 
collapse in regions where they are magnetically supercriti- 
cal, while their subcritical fragments may remain supported 
for times significantly longer than their free-fall time. 

One fundamental distinction between the hypotheses 
of local turbulence and of LSI is that, in the former, 
the kinetic energy of the turbulent motions internal to a 
clump is assumed to act fully as support against gravity, 
while in the latter, part of the kinetic energy of these 
motions may be globally compress ive, either promoting 



collapse or being a consequence of it (Hunter fc Fleckl 
Ballcstcros-Parcdes, Vazqucz-Scmadcni, fc Scald 



1982 



199c 



Ballesteros-ParedesI [2OO6I: iDib et aIll2007^■ Recent models 



of th e SFR (e.g., lElmegreenI |2002| : iKrumholz fc McKed 
120051 ') or of the turbulent clump mass f unction as the origin 
of the stellar initial mass functio n (|Padoan fc Nordlundl 
l2002l : iHennebelle fc Chabrieij I2OO8I ) have been formulated 
under the assumption of local turbulence. It is thus impor- 
tant to determine whether this assumption is verified in 
numerical simulations of self-gravitating turbulence. The 
realization of the LSI scenario would thus imply a reduced 
amount of support against self-gravity. 

Moreover, t he models by Padoan fc Nordlundl l|200a, 
hereaf ter PN02'1. lElmegreenI (|2002l ) and lKrumholz fc McKcc' 



20051 . h ereafter KM 5) rely on the idea advanced byi Padoan 
19951), lElmegreenI (|2002l ') and IVazguez-Semadeni et al.l 



2003) that the mass that proceeds to collapse is that 



which is deposited by the turbulence in subsonic, yet Jeans- 
unstable ( "super- Jeans" ) fragments. The latter authors pro- 
vided indirect evidence that this could be so by showing the 
existence of a correlation between the SFE and the so called 
"sonic scale" of the turbulence, the scale below which the 
turbulent motions are subsonic on average. However, it is not 
necessary that only the mass in these subsonic, super- Jeans 
structures proceeds to collapse. Material in supersonic, yet 
effectively gravitationally unstable structures can also par- 
ticipate in the collapse. 

The goal of the present paper is to contribute towards 
the understanding of the role of the velocity field's topology 
on the control of the star formation process. To this end, 
we use numerical simulations (described in!}21) of randomly 
driven turbulence in isothermal, self-gravitating flows aimed 
at investigating whether the hypotheses and predictions of 
analytical models are verified in the non-magnetic case. The 
simulations have different sizes, mean densities and veloc- 
ity dispersions, but scaled so that all have approximately 
the same value of the virial parameter a = 2_Ekin/|£^grav|, 
thus satisfying the hypotheses of the KM05 model. In these 
simulations, we first search for the fraction of super- Jeans, 
subsonic subregions in a cloud, to see whether they can be 
deemed responsible for the mass that ends up collapsing 
(33I- Second, we investigate whether the internal velocity 
dispersion in dense subregions of a cloud can be considered 
as random as that at larger scales O, or instead exhibits 
increasing amounts of the compressive component as the 
density of the structures increases and their size decreases. 
In !j5l we ask whether the suite of simulations agrees with 



the prediction of the model by KM05 for the dependence 
of the SFRff on the turbulent Mach number. Finally, in !j6] 
we summarize and discuss the implications of our results, in 
particular comparing with previous work. 



2 THE MODELS 

We have performed three simulations of non-magnetic, self- 
gravitating, isothermal turbulence at a resolution of 512^ 
zone s, using a to tal variation diminishing (TVD) scheme 
( Kim et al.l 1 19991 ) with periodic boundaries and random 
Fourier driving with a spectrum P{k) = fc^ exp(— 8fc/fcpk). 
Here, fc is the wavenumber and fepk = 2(27r/L) is the energy- 
injection wavenumber, with L being the computational box 
size. The energy is thus injected mostly at scales of or- 
der half the box size. The driving is purely rotational (or 
"solenoidal" ) , thus having no compressive component. A 
prescribed rate of energy injection is applied in order to 
approximately maintain the rms Mach number Ms = cr/cs 
near a "nominal" value that characterizes the run. Here, a is 
the three-dimensional velocity dispersion and Cs is the sound 
speed, taken equal to 0.2 km s~^ in all runs (corresponding 
to T = 11.4 K). The actual value of Ms fluctuates and is 
slightly different from the nominal value, because the nu- 
merical scheme is designed to maintain a constant energy 
injection rate, not a constant rms Mach number. 

The other parameter that characterizes a simulation is 
the "Jeans number" J = L/Lj, where Lj = {tvc^/GpY'^ is 
the Jeans length, with p = pno being the mean density of 
the simulation, no the mean number density, p. = 2.36 thh 
the mean particle mass, and mn the mass of the Hydrogen 
atom. All three simulations are evolved for 3.2 turbulent 
crossing times before turning on self-gravity, in order for 
the turbulence to reach a fully developed state, and thus 
avoid applying the self-gravity directly on the imprints of the 
random driving. The three simulations differ in physical size 
L, mean density no and nominal rms Mach number Ms, but 
in each case their values are chosen as to give the same value 
of the ratio Ms/ J. The nominal values of the pairs (Ms, J) 
for the three runs are (8, 2), (16, 4), and (24, 6), respectively, 
corresponding to a nominal value of Ms/J = 4. The runs 
are named mnemonically by means of their values of Ms and 
J. 

Note that {Ms/ J)'^ is proportional to the virial param- 
eter a, as can be seen by approximating -Bkin ~ Mcr'^ /2, 
where M is the total mass in the simulation, and |-Egrav| ~ 
GM^ / L. For a spherical cloud of mass M = -iTrpL^ /3, we 
thus have 



2-Ekin 

-L/crraV 



GMyi 



3 Ms^ 
4^72"' 



(1) 



Our simulations thus all have a nominal value of q ~ 1.22. 
The actual average values of Ms and of a, together with 
other parameters of the runs, are indicated in Table [2] in 
which the column labeled tgrav gives the time at which self- 
gravity is turned on, and the column labeled "Ms (real)" 
gives the actual measured average value of the rms Mach 
number over the duration of the run. 

It is also worth noting that fixing the ratio Ms/J only 
fixes the ratio a /{p^'^L) (at a given Cs), and so we still have 
freedom to choose the values of the individual physical pa- 
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Run parameters 



Name 


L 


"0 


M 


i,7 


J 


Ms (real) 


a (real) 


"^rms 


<ff 


igrav 


tturb (nominal) 


grid cell size 




[pc] 


[cm-3] 


[Mq] 


[pc] 








[ km s-i] 


[Myr] 


[Myr] 


[Myr] 




[pc] 


Ms8J2 


1 


2000 


115.8 


0.5 


2 


8.6 


1.4 


1.7 


2.5 


2 


0.625 




0.00195 


Msl6J4 


4 


500 


1853 


1 


4 


15.7 


1.2 


3.1 


5 


4 


1.25 




0.00781 


Ms24J6 


9 


222.22 


9382 


1.5 


6 


23.0 


1.1 


4.6 


7.5 


6 


1.875 




0.0175 



rameters. We do this by assuming that the size L and mean 
density no of our simulations satisfy one of Larson's (1981) 
relations, i.e., no oc L~^ . Since all three simulations have 
the same nominal value of a, then the above assumption 
also implies that our suite of simulations also satisfies the 
other Larson relation, a oc L^'^. The set of physical values 
for the simulations are also reported in Table (2] 



3 FRACTION OF SUBSONIC, SUPER-JEANS 
STRUCTURES 

3.1 Procedure 

In this section, we measure, as a function of region size, the 
fraction of regions in the numerical simulations that is both 
subsonic and super- Jeans, in order to test the hypothesis 
that these are indeed the structures that collapse gravita- 
tionally to form stars. For generality, we consider two types 
of regions: a) the set of all cubic sub-boxes of a simulation 
of a given size and b) dense clumps defined by a density 
threshold criterion. The first set contains both overdense 
and underdense regions, and is intended to provide unbiased 
statistics of the velocity field in all subregions of a given size 
within the numerical box, independently of the local density. 
The second set contains only overdense clumps. 

To isolate the efi'ect of self-gravity, we perform the pro- 
cedure at two different times in each simulation. First, we 
consider the last data dump before gravity is turned on, at 
which the density distribution must be a consequence of the 
turbulent flow alone. Second, we consider a data dump at 
around two free-fall times ts after having turned gravity on, 
at which significant gravitational collapse has occurred, and 
the density structure should be the result of the combined 
effects of turbulence and self-gravity. 

For the analysis using sub-boxes of the simulation, we 
subdivide the latter in cubic sub-boxes of sizes 2, 4, 8, 16, 32, 
64, and 128 grid zones per dimension. Since the simulations 
are performed at a resolution of 512 grid cells per dimension, 
and runs Ms8J2, Msl6J4 and Ms24J6 respectively have sizes 
L = 1, 4 and 9 pc, the grid cell size differs for each run, and 
is also indicated in Table [2] It should be noted that, since 
the sub-boxes are located at fixed positions within the simu- 
lation box, their locations bear no special relation with those 
of actual clumps. A clump may be located at the interface 
between two sub-boxes, and it in general moves from one 
sub-box to another as time elapses. However, we make no 
attempt to repeat the procedure shifting the positions of the 
sub-boxes to find a possible "best" match between the sub- 



boxes and the actual clumps. We consider that the number 
of sub-boxes is large enough that no particular choice of the 
origins of the boxes is better than any other. 

For the analysis with clumps, we define a clump as a 
connected set of grid cells with densities above a certain 
threshold. To create an ensemble of clumps, we consider a 
series of thresholds, of n = 32, 64, 128 and 256 times the 
mean density of the simulation. We then approximate their 
size by L = (31^/47r)^'^, where V i s their volume. As dis- 
cussed bv IVazguez-Semadeni et al.l l|2005ah . this is a rather 
robust estimator of the clumps' size. Finally, for the plots, 
the clumps are classified by their sizes (independently of the 
threshold from which they originated) and added to loga- 
rithmic size bins defined by successive powers of 2 from 2 to 
64 grid cells. 

For each region (sub-box or clump), we compute its 
mean density, its internal three-dimensional velocity dis- 
persion u^ (substracting the region's mean velocity), and 
whether it is Jeans stable or unstable. In the case of sub- 
boxes, since their size is predetermined and the tempera- 
ture is constant in space and time, we simply compute the 
critical density for Jeans instability (i.e., with respect to 
thermal support onlyjj at the specified size. Since the flow 
is isothermal, this "Jeans density" is unique at a given scale, 
and thus we can compare it with the sub-box's mean density. 
For the clumps, we simply compare their estimated size to 
the Jeans length derived from the clump's mean density. Fi- 
nally, we count the number of regions that have a size larger 
than their associated Jeans length, the number of regions 
that have a subsonic velocity dispersion, and the number of 
regions that satisfy both conditions simultaneously. 

Note also that, in the case of clumps, we have not cor- 
rected for the possibility of them being located at the bound- 
aries of the periodic numerical box. So, if a clump crosses 
one or more boundaries, it is artiflcially split into two or 
more fragments by our clump-defining algorithm. However, 
we do not expect this minor omission to introduce any seri- 
ous biases since, as we discuss in i]3.2l the fraction of super- 
Jeans clumps is found to already be unity at the largest 
sizes (clumps of even larger size would be even more super- 
Jeans) , and the fraction of subsonic clumps is already zero 
at scales significantly smaller than the largest clump sizes 



We do not consider the non-thermal contribution to the "sup- 
port" , since we are only interested in the fraction of simultane- 
ously subsonic and super-Jeans structures. Moreover, as discussed 
in i|4] it is not clear that all of the non-thermal kinetic energy can 
be considered to provide support. 
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Figure 1. Fraction of subsonic [triangles, dotted lines) and of super-Jeans (diamonds, solid lines) sub-boxes in runs Ms8J2 (top row), 
Msl6J4 (middle row), and Ms24J6 (bottom row), as a function of sub-box size. The entire numerical box is subdivided in sub-boxes of 
the indicated size. The left panels show the fractions shortly before the time when self-gravity is turned on (tgrav). The right panels show 
the fractions at approximately two free-fall times after igrav • The fraction of sub-boxes that are both subsonic and super- Jeans is zero 
at all sub-box sizes, and thus cannot be shown in this figure. 



found (clumps of larger size would be expected to be even 
more supersonic). Besides, the sizes of the clumps are never 
larger than l/lOth of the numerical box, and so the proba- 
bility of them crossing the boundary is relatively low. 

3.2 Results 

The fractions of subsonic and super- Jeans structures for the 
three runs are shown for sub-boxes in Fig. [l] and for clumps 



in Fig. [21 In both figures, the top row shows results for run 
Ms8J2, the middle row shows run Msl6J4, and the bottom 
row shows run Ms24J6; the left panels show the result before 
the time at which self-gravity is turned on (igrav; cf. Table[5)), 
and the right panels show it at approximately two free-fall 
times after igrav Note that, for the sub-boxes, the fractions 
can be very small, and are thus shown in logarithmic scale. 

The most notable result of this analysis is that the set 
of simultaneaously subsonic and super- Jeans structures (ei- 
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Figure 2. Same as Figure [T]but for clumps rather than sub-boxes. The clumps are defined as connected structures with densities above 
a density threshold. The ensemble of clumps was created by considering thresholds 32, 64, 128 and 256 times the mean density no. The 
fraction of clumps that are both subsonic and super- Jeans is again zero at all clump sizes considered (not shown). Scales at which there 
are no clumps have no points drawn. 



ther sub-boxes or clumps) is empty at all the scale sizes we 
sampled. Thus, this fraction is not plotted in Figures [l] and 
[21 We discuss the implications and limitations of this result 
in iH 

Some other features are worth noting. The fraction of 
subsonic {triangles, dotted lines) sub-boxes or clumps as a 
function of size shows no clear trend with the inclusion of 
self-gravity. However, the fraction of super- Jeans {diamonds, 
solid lines) structures at small (for sub-boxes) or intermedi- 
ate (for clumps) scales tends to increase in the presence of 



self-gravity. In fact, for the small-scale run Ms8J2 there are 
no super- Jeans sub-boxes in the absence of self-gravity, but 
significant amounts appear after it has been turned on. This 
means that the presence of self-gravity changes the distri- 
bution of sub-box masses in comparison to that produced by 
turbulence alone, increasing the fraction of regions that can 
proceed to gravitational collapse. That is, the effect of self- 
gravity can begin prior to the actual "capture" of a region to 
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proceed to collapse|j In other words, a density enhancement 
is in general expected to reach higher peak and mean densi- 
ties in the presence of self-gravity than when it is produced 
only by turbulence. We refer to this as gravity aiding in the 
production of the density fluctuation, regardless of whether 
it eventually collapses or not. 

This conclusion is also supported by the probability 
density function (PDF) for the density field of the three 
runs. In Fig. [3] we show the density PDFs of the three runs 
at 0, 1 and 2 times the global free-fall time, tg, of each sim- 
ulation (cf. Table [2)). The PDF at tg is representative of 
the effects of turbulence alone, while those at 1 and 2 iff 
show the efi'ect of turbulence and gravity combined. In all 
runs, the PDFs in the presence of self-gravity show a promi- 
nent high-density tail, implying that the relative frequency 
of high density regions is higher in this case, compared to the 
effect of turbulen ce alone. Similar results have been found 
by other workers (|Klessenll2000l : lDib fc BurkertlbOOSl ). Thus, 
the production of super- Jeans structures is itself aided by the 
inclusion of self-gravity. 

One final point to note is that, since higher density 
thresholds in general imply that the resulting clumps have 
smaller sizes and higher densities, the small clumps are gen- 
erally also denser. Thus, Fig. [5] implies that there indeed 
exists a population of dense, subsonic clumps in the simula- 
tions, which, however, are not superjeans. 



4 COMPRESSIVE COMPONENT OF THE 

VELOCITY FIELD IN DENSE STRUCTURES 

We now investigate the nature of the velocity field in subre- 
gions of our turbulent supersonic flows, in particular aiming 
at whether our small-scale simulation Ms8J2 is statistically 
representative of regions of similar size within the large-scale 
run Ms24J6. For this purpose, we again subdivide the large- 
scale simulation's (Ms24J6) domain into sub-boxes, but this 
time having the same size as the small-scale run Ms8J2. 
Note, however, that run Ms8J2 has a size of l/9th that of run 
Ms24J6. Since the latter has a resolution of 512 grid cells per 
dimension, then the sub-boxes should span 512/9 = 56.89 
cells per dimension, which we round to 57 cells. Due to this 
rounding, we cannot fit 9 full, non-intersecting sub-boxes 
along each direction within the whole box, and thus we only 
consider the first 8 sub-boxes from the origin in each direc- 
tion, for a total of 8^ = 512 sub-boxes, each with a linear 
size of 57 grid cells. 

For each sub-box, we measure its mean density and the 
mean divergence of its velocity field. To compute the diver- 
gence, we Fourier-transform the three velocity components 
for the entire simulation box, and compute the divergence 
in Fourier space as 



J^CV ■ v) = —ik ■ Vk 



(2) 

where i — \/—l, J-{) denotes the Fourier transform of its 
argument, k is the wavevector, and Vk is its associated 

^ We say that a certain density enhancement is "captured" by 
gravity when it was originally produced by a turbulent compres- 
sion, but at some time while its mass and/or density are increas- 
ing, it suddenly becomes gravitationally unstable and begins to 
undergo gravitational collapse. 



Fourier velocity amplitude. We then transform back to phys- 
ical space to obtain the divergence field for the whole box, 
and finally we take the average of this field in each sub-box. 

Figure |4] shows the result of this exercise, giving the 
mean divergence of each sub-box of run Ms24 J6 as a function 
of its mean density at two times before self-gravity is turned 
on: at f = 4 Myr {— 2.1 fturb, where tturb = L/vrms ~ 1.875 
Myr is the "nominal" turbulent crossing time; cf. Table [5} 
on the left panel, and at t = 6 Myr = 3.2 tturb on the right 
panel. The contours show the two-dimensional histogram 
of the points, and are drawn at levels 0.143, 0.286, 0.428, 
0.571, 0.714 and 0.857 of the maximum. Figure [5] shows the 
corresponding plot at times in which self-gravity has been 
on for 1 iff {left panel) and 1.93 ta {right panel), where ta = 
Lj/cs is the free-fall time. The vertical lines in both figures 
show the mean densities of runs Ms24J6 {dotted line) and 
Ms8J2 {dashed Ime). 

In all cases, the contours are seen to have an elongated 
shape, indicating that, although with abundant scatter, a 
positive correlation exists between the mean value of —V ■ v 
(i.e., the velocity convergence) and the mean density of the 
sub-boxes. The contours are more elongated in the cases 
with self-gravity, for which we find an average fit of 



V- V 



km s ^pc ^ 



-(0.45 ±0.05) logjo 
0.13 ±0.03, 



222 cm-3 



(3) 



where the uncertainties in the slope and the intercept span 
the range of values we have encountered in the two cases 
sampled. The fit implies that, at the size and mean density 
of run Ms8J2, which is 9 times smaller and denser than run 
Ms24J6, a mean divergence of ~ —0.6 km s~^pc~^ should 
be expected. Instead, run Ms8J2, which is driven by random 
forcing at its own scale in order to mimic the assumption of 
local turbulence, has zero mean divergence by construction. 
The trend of V • u with mean density observed in the sub- 
boxes of the large-scale simulation Ms24J6 suggests that the 
LSI scenario (cf. 2} *■' verified in the density enhancements 
even m driven turbulence regimes. 

It is also noteworthy that the probability of producing 
sub-boxes of the same density and size as run Ms8J2 in run 
Ms24J6 is extremely low (1 case in 1024, since each plot 
contains 512 points) in the absence of self-gravity (Fig. |4]). 
This probability increases substantially in the presence of 
self-gravity (13 cases in 1024). Thus, we conclude that the 
formation of regions with the same size and density as those 
of run Ms8J2 within run Ms24J6 seems to require the pres- 
ence of self-gravity. This result further reinforces the con- 
clusion from |3]that self-gravity intervenes in the formation 
of super- Jeans structures, and not only in their collapse. 

It is interesting to compare the magnitude of the con- 
vergence seen in Figs. U and [S] with that of the vorticity, in 
order to get an idea of which type of motion (potential or 
solenoidal) dominates within the subregions of run Ms24J6. 
Again, we calculate the vorticity in Fourier space as 



T{\/ X v) = —ik X Vk, 



(4) 



and use a similar averaging procedure for the Ms8J2-sized 
sub-boxes of run Ms24J6 as for the divergence. Since the 
vorticity is a (pseudo-)vector, in Fig. [6] we plot the magni- 
tude of the average vorticity in each sub-box. This is shown 
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Figure 4. Negative mean divergence of the velocity field as a function of mean density for sub-boxes of the large-scale simulation Ms24J6 
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Ms24J6. 



at i = 4 Myr (a time without self-gravity; left panel) and at 
t = 14 Myr (a time with self-gravity; right panel). In gen- 
eral it is seen that the magnitudes of the divergence and of 
the vorticity are comparable, suggesting that the sub-boxes 
contain comparable amounts of the two types of motion. 

Note, however, that it is not possible from these figures 
to obtain a quantitative estimate of the amount of kinetic 
energy contained in each type of motion, since, contrary to 
the case of the energies, not the divergence nor the vortic- 
ity nor their sum are bounded or conserved in any sense. 
It would be necessary to compare the kinetic energy in each 
type of motion directly but, because this task is not straight- 
forward, we defer it to a future study. In any case, our 
simulations show that, to order of magnitude, the potential 
and solenoidal parts of the velocity gradient are comparable. 



and therefore we expect the fraction of compressive energy 
present in dense regions to be substantial. 



5 STAR FORMATION EFFICIENCY IN 
CONSTANT- VIRIAL-PARAMETER 
STRUCTURES 

In this section we now proceed to measure the SFE in 
terms of what KM05 called the star formation rate per free- 
fall time, SFRff, in the three simulations, in order to test 
whether the predictions of their model are verified. Follow- 
ing KM05, we define the SFRff as the fraction of a cloud's 
mass that is deposited int o stars in one free-fall time. As 
in previous papers (e.g.. IVazguez-Semadeni et aU l2005bl : 
iGalvan-Madrid et al.ll2007h . we approximate this fraction 
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Figure 5. Same as Fig. U] but at two times with self-gravity: i = 14 Myr (Je/t panel) and t = 20 Myr {right panel) (respectively 1 iff 
and 1.93 tg after having turned on self-gravity, with iff being the global free-fall time of the simulation [run Ms24J6; see Table[2]). The 
slopes of the fits through the data points have slopes 0.50 {left panel) and 0.40 {right panel). We note that in the presence of self-gravity, 
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latter is significantly increased with respect to the case without self-gravity (Fig.l4t. 
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Figure 6. Same as Figs. |4] and [5] but for the magnitude of the vorticity at t = 4 Myr {left panel, without self-gravity) and i = 14 Myr 
{right panel, with self-gravity) (respectively, compare to the left panels of Figs. |4] and |5]l. It is seen that comparable magnitudes of the 
(negative) divergence and of the vorticity occur before and after turning self-gravity on. 



as the fraction of the total mass in the simulation that 
is in regions with density above a certain threshold. The 
threshold is chosen to be unambiguously indicative that col- 
lapse has occurred. Specifically, we take a threshold density 
nthr = 1000 no, which is a much larger density than can 
be achieved through the turbulent fluctuations alone in any 
of the runs, and guarantees that the mass contained in 
objects above those densities is either collapsing or gravita- 
tionally locked within them and does not redisperse (except 
for sligh t fluctuations, discussed in the nex t paragraph jj. As 
noted in lVazguez-Semadeni et al.l (|2005a ), the collapsed (or 
"accreted") mass fraction does not depend heavily on the 
chosen threshold, as long as the latter satisfies the above 



^ Note that, in our simulations, gravitational collapse is necessar- 
ily terminated when most of the collapsed mass is deposited into 
a single or a few grid cells. It is in this sense that the material 
is "locked" into the collapsed objects, and by no means it should 
be interpreted as implying that a real hydrostatic object (other 
than one or several stars) has formed. 



conditions. As a confirmation, we have also used a thresh- 
old nthr ~ 500 no, and found no significant difference in the 
results, as described below. Note however that our procedure 
may overestimate the actual mass deposited into protostel- 
lar objects, since it assumes a 100% efficiency of conversion 
of gas into stars above the threshold density, while in ac- 
tual cluste r-forming cores, the expected efficiency is closer 
to 30-50% (|Lada fc Ladall2003h . 

The left panels of Figure [7] show the evolution of the 
mass accreted onto collapsed objects as a function of time 
after having turned gravity on, in units of Myr {upper frame) 
as well as in units of the free-fall time iff {lower frame). The 
SFRff is then the slope of the collapsed mass fraction versus 
time, the latter in units of iff. 

The accretion histories of the runs are seen to be noisy, 
even exhibiting slight, occasional drops. This refiects the 
fact that our numerical scheme does not include any pre- 
scription for sink particles or cells, and therefore the mass 
in the collapsed object is not fully "locked" in the object - 
the outer parts of it may have densities that oscillate about 
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Figure 7. Le/i panels: Fraction of mass accreted onto collapsed objects as a function of time for the three simulations, with time in 
units of Myr {top left) and of the free-fall time of each simulation {bottom left). The average SFRji is computed by fitting a least-squares 
line to the plot of accreted mass vs. time-in- units-of-iff. Right panel: Star formation rate per free-fall time (SFRff, i.e., average accreted 
mass after one free-fall time) as a function of the time-averaged rms Mach number of the simulation for the three runs. Triangles (resp. 
squares) denote results obtained with a density threshold n^hr = 1000 no (resp. ntjjj = 500 no) for defining the collapsed objects. The 
dotted and solid lines show the best fits through the data points for nti^^ = 1000 no and n^jir = 500 no, respectively. The slope of the 
dotted line is —0.55, while that of the solid line is —0.58. The dashed line shows the prediction of the KM05 model, with slope —0.32. 



the threshold we use for defining the collapsed objects. This 
implies that the SFRh is determined only to within a certain 
uncertainty in our simulations. For estimating the average 
SFRff at each Ms , we thus fit a least-squares line to the plot 
of accreted mass vs. time (the latter in units of iff; lower 
left panel of Fig. O, and show error bars corresponding to 
±3(7, where a is the standard deviation of the fit. Note also 
that the accreted mass for run Ms8J2 appears to saturate at 
t ~ 5 Myr ~ 2ts, so for this run we only take into account 
times shorter than this saturation time. 

With this procedure, we obtain the average values of the 
SFRff shown in the right panel of Figure [T] For reference, 
we show the results obtained with the two thresholds nthr = 
500 no and nthr = 1000 no. We see that the data points 
are shifted slightly to higher values for the lower threshold, 
since the structures contain more mass, but the trend of the 
SFRff is essentially the same in both cases. 

Our results can be compared with the prediction of the 
KM05 model, which those authors fitted by 



SFRff ^ 0.014 



(n) 



/M^X-O-32 

ViooJ 



(5) 



Thus, a residual dependence on the Mach number ~ M^ 
is expected in the KM05 model even at constant a. The slope 
of —0.32 is indicated in the right panel of Figure [7] by the 
dashed line, while the solid and dotted lines indicate fits to 
our data, with slopes —0.58 for nthr = 500 no and —0.55 for 
nthr = 1000 no, respectively. 

We see that, within the uncertainties, our results can ac- 
comodate the slope predicted by KM05, although the agree- 
ment is only marginal, because the most probable slope that 
we measure is significantly steeper than theirs. Besides, the 



intercept of the fitted line appears to imply a larger coef- 
ficient than the one in eq. (O. For example, at Ms = 10, 
the dotted line (nthr = 1000 no) in the right panel of Fig. 
Ogives SFRft~ 0.15, while eq. ([5]), with a = 1.2, predicts 
SFRff = 0.031. 

Several effects, both physical and of procedure, can ac- 
count for these differences. Possible effects of procedure are: 
a) We are assuming a 100% efficiency of conversion from 
gas to stars above the density threshold. Instead, this ef- 
ficiency is known to be ^ 3 0% for high-mass star form- 
ing cores (|Lada fc Ladall2003i ) and even lower for low-mass 
ones, so we are clearly overestimating the total mass de- 
posited into stars by a factor of 3-10|j b) KM05 calibrated 
their model using linear appr oximations to the SFEs of 
IVazguez-Semadeni et al.l (|2003l ). while in reality the mass 
accretion histories presented by those authors are strongly 
nonlinear, c) We have defined the free-fall time simply as 
iff = I/j/cs, which is larger by a factor of •y/32/37r than 
the alternate definition tjf = (37r/32Gp)^". Our choice is 
a reasonable estimate, since actual collapse appears to take 



* Note that, because we use a threshold that is simply a fixed 
multiple of the simulation's mean density, we may be biasing the 
measured slope, since the physical threshold density is highest 
in the small-scale run Ms8J2 and lowest in the large-scale one, 
Ms24J6. However, correcting for this by, say, taking a fixed phys- 
ical density in all three simulations will actually tend to make 
the slope of the SFRff vs. M^ curve even steeper, since a higher 
threshold in code units should be used in Ms24J6, leading to an 
even smaller measured SFRff. Thus, this effect cannot account 
for the fact that we measure a steeper slope than that predicted 
by KM05. 
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from ~ 1.6 to a few times t{f ( LarsonI 1 19691 : iGomez et al.1 
I2OO7I : iGalvan-Madrid et al.ll2007r ). due to the non-negligible 
role of the thermal p ressure. However, KM05 used tg 
l|Krumholz fc Tarj|2007l ). so our measured masses are larger 
than those that would be measured out to iff by a factor of 
a few, raising the intercept of our fit with respect to that of 
KM05. 

On the other hand, possible physical efi'ects are: a) the 
fact discussed in |4]that the clumps in each simulation lack 
support in comparison with the KM05 model because the 
motions are not fully random, but instead involve a signifi- 
cant inflow component, b) Most importantly, the mass that 
appears to be collapsing in our simulations is not that in 
subsonic, super- Jeans regions, but rather the mass in re- 
gions in regions that are supersonic but still gravitationally 
unstable. This allows for a larger amount of mass involved 
in the collapse, since in general regions with larger velocities 
are expected to be larger, and thus more massive. In conclu- 
sion, it is not possible to determine whether the differences 
between the KM05 prediction for SFRff and our results are 
significant. However, it is clear that the mechanism of col- 
lapse a cting in o ur simulati ons is different from what t hey, a s 
well as IPadoanI (119951 ') and IVazguez-Semadeni et al.l (|2003n , 
assumed, since there are no subsonic, super- Jeans structures 
in our simulations. 



6 DISCUSSION AND COMPARISON TO 
PREVIOUS WORK 

6.1 Velocity convergence 

Our results have a number of implications for our under- 
standing of the velocity field in MCs and for analytical mod- 
els of star formation. A first result is that in sub-boxes of the 
large-scale run (Ms24J6) of the same size as the small-scale 
run (Ms8J2), the velocity field exhibits a clear trend towards 
being convergent if the subregion is overdense with respect 
to the rest of the simulation. This suggests that run Ms8J2, 
with its zero overall mean divergence, is not representative 
of the typical region of the same size embedded within a 
larger medium, since it lacks an expected mean divergence 
of ~ —0.6 km s^^pc"^, as predicted by eq. Q. 

This result also suggests that a significant compo- 
nent of the observed linewidths in MCs and their sub- 
structure should be compressive (i.e., conform with the 
LSI scenario), in agreeme nt with the original suggestion of 
iGoldreich fc KwanI ([197J) . It is even more relevant that this 
result is observed in numerical simulations of driven turbu- 
lence (with the driving applied at the large scales and being 
purely rotational), contrary to the frequent belief that such 
large-scale inflows occur only in decaying turbulence. It is 
necessary to stress, however that our analyses cannot dis- 
tinguish whether the convergent motions in the overdense 
regions are a cause or a consequence of gravitational con- 
traction. This is beyond the diagnosing capabilities of the 
studies we have performed. 

Dimensionally, the values of the mean divergences we 
have found for the sub-boxes of our large-scale simulation 
can be comp ared to the typical valu es of the velocity gradi- 
ents found bv lGoodman et aU (| 19931 ) in dense cores of MCs, 
which they however interpreted as representative of uniform 



rotation. For cores with typical sizes between ~ 0.1 and 1 pc, 
those authors found velocity gradients ranging between 0.3 
and 4 km s~^ pc~^. Our average value of the velocity diver- 
gence for regions of the same size and density as those of run 
Ms8J2, V • t; ~ —0.6 km s~^ pc~^, seems to be on the low 
side of this distribution, although it may seem reasonable if 
we consider that only a fraction of the kinetic energy should 
be compressive (i.e., non-rotational) in general. We can also 
compare with Larson's (1981) linewidth-size relation, which, 
for ^^ CO and "CO data reads Av » 1 km s-'^ \L/l pcj^/^ 
(e.g.. ISolomon et al.|[l987l : iHever fc Brunt! |2004| ). Thus, for 
L ~ 1 pc, Larson's relation predicts a velocity dispersion 
roughly 1.8 times that of our "typical" velocity conver- 
gence at the scale of run Ms8J2, given by eq. ((3}, suggesting 
that, on average, the velocity dispersion contains compa- 
rable amounts of compressive and non-compressive energy. 
Of course, the fluctuations are large on both our velocity 
convergence-density relation (eq. |3} and on Larson's rela- 
tion, so large deviations from this mean trend can be ex- 
pected in individual clumps. 

Since the compressive part of the velocity works to 
promote compression, the main implication of our result 
is that not all of the non-thermal kinetic energy in clouds 
and clumps is available for support against gravity, a fact 
that has been overlooked by analytic models of star forma- 
tion from the turbulent conditions in molecular clouds and 
clumps. In particular, the models by PN02 and KM05 as- 
sume that a star-forming core must exceed t he thermal Jeans 
mass in order to collapse but, as shown bv lHunter fc Fleckl 
( 1982), this mass is reduced i n the presence of a compressive 
velocity field. The model by iHennebefle fc Chabrierl (|2008l ) 
does consider an additional "support" by the turbulence, 
but this may again be reduced if the "turbulence" actually 
contains a significant fraction of inwards motions. 

The LSI scenario is naturally motivated by the 
fact that only compressive motions can produce density 
fluctuations, while vortical modes are incompressible. 
Thus, if the density enhancements (clumps) in a flow 
are produced by turbulent fluctuations, the velocity field 
within these clumps is likely to still exhibit the signature 
of the external convergent motions that form ed them 
(|Ballesteros-Paredes. Vazguez-Semadeni. fc Scald |l999J). 
This scenario also has the implication that the largest veloc- 
ity dispersions are expected to occur not in the densest parts 
of the structures, but at their outskirts, sin ce the densest 
gas h a s been shocked and has slowed down (JKlessen et al.l 
I2OO5I : IComez et all |2007| ). This effect has been found 
observationally, albeit with lower velocity dispersion than 
in the models (see, e.g.. lAndre. Basu fc Inutsukall2008l . and 
references therein). Finally, the LS I scenario is also fully 
consistent with the observation by iHever fc BruntI (120071 ) 
that Principal Component Analysis of the velocity field 
in molecular clouds and their substructure systematically 
shows the dominance of a whole-cloud dipolar pattern. 

At first sight, one could think of two possible ways 
to avoid the LSI picture. One possibility would occur if 
the clumps were quasi-static struct ures in single-phase me- 
dia confined by ram pressure ([Bertoldi fc McKee 1993). 
In this case, the boundaries must be accretion shocks 
(JFolini fc WaldeJl2006l : IWhitworth et al.ll2007l : IComez et all 
200^), so that the mass of the structures must grow over 
time, possibly eventually becoming strongly self-gravitating 
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and proceeding to collapse. Thus, even though such clumps 
can be quasi-static in their central parts, they must be sur- 
rounded by an accreting envelo pe that involves a net conver- 
gence of the velocity field f^Gomez et al.1 120071 ). The other 
possibility would occur if the med ium, even within MCs, 
turns out to be thermally bistable (JHennebelle fc Inutsukal 
120061 ). in which case the clumps could be hydrostatic 
objects bound by the thermal pressure of their warm, 
te nuous environment, as in th e class ic two-phase model 
of iField. Goldsmith, fc Habind (|l969l ). However, even if 
molecular clouds are thermally bistable, clumps that are 
at much higher thermal pressures than their surroundings 
must be driv en either by ram-pre ssure compressions or by 
self-gravity (|Baneriee et al.l 120081 ). So, there appears to be 
no escape from the need to have convergent flows involved 
in the formation and evolution of the densest, star-forming 
clumps. 



6.2 Absence of subsonic, super- Jeans structures 

Another result we have obtained is that, in our driven- 
turbulence simulations, no simultaneously subsonic and 
super- Jeans structures were found, either before or after self- 
gravity was turned on. Of course, this result does not rule 
out the existence of such structures, and in f act subsonic , 
super- Jeans cores are routinely observed (e.g., iMverdl 19831 : 
I Andre et al.ll2007l ). Our result may be an artifact of the tur- 
bulence being continually driven and/or the absence of mag- 
netic fields and/or the fact that the actual values of a in our 
simulations are somewhat greater than unity. Nevertheless, 
our simulations produced abundant collapse, indicating that 
the format ion of simultaneously subsonic and super- Jea ns 
structures (|Padoanlll995l : IVazguez-Semadeni et al.ll2003l ) is 
not the only possible route to collapse. Since this notion is at 
the foundation of models such as those of PN02 and KM05, 
it is likely that those models may need to be revised to con- 
sider the possibility that stars may form via the collapse 
of larger-scale, supersonic regions, in which the motions are 
not fully supportive against gravity. 

It is also importa n t to reconsider the results of 
IVazguez-Semadeni et al.l (120031 ) in the light of our present 
results. In that paper it was shown that there exists a cor- 
relation between the sonic scale of the turbulence and the 
SFE, so that, as the sonic scale becomes smaller (at con- 
stant J), the SFE decreases. This was interpreted as indi- 
rect evidence that the available mass for collapse decreased 
as the sonic scale became smaller, and therefore that the col- 
lapsed objects indeed originated from subsonic, super- Jeans 
structures. However, our finding that there exist alterna- 
tive routes to collapse that are not based on the formation 
of subsonic, super - Jeans structures suggest s that the cor- 
relation found by IVazguez-Semadeni et al.l (|2003l ) may be 
simply indicative of a general scaling of both the SFE and 
the sonic scale with rms Mach number, but not that the 
fraction of mass in subsonic, super- Jeans structures directly 
measures the mass that is on route to collapse at any given 
time. 



6.3 Role of gravity in the formation of 

Jeans-unstable structures and the density 
PDF 

Our simulations also suggest that self-gravity is not only 
involved in the capture of turbulent density fluctuations to 
make them collapse, but also in the production of collaps- 
ing objects, as shown by the distortion of the density PDF 
and by the increase of the fraction of self-gravitating regions 
at small scales in the presence of self- gravity. In this sense, 
these results may repre sent the driven-t u rbulen ce counter- 
part of the results by IClark fc Bonnelll l|2005l ). who con- 
cluded from decaying turbulence simulations that the tur- 
bulence does not directly produce the collapse of clumps by 
making them reach their own Jeans mass, but rather just 
produces the seeds for subsequent gravitational fragmenta- 
tion of the large-scale gravitationally unstable structures. In 
our driven case, it appears that turbulence alone produces 
only a few super- Jeans structures, while, in the presence of 
self-gravity, super- Jeans objects are much more readily pro- 
duced. 

This result is likely to have an implication for the PN02 
model. This model assumes that turbulence alone is re- 
sponsible for the formation of the cores, with gravity only 
playing a role if they become Jeans unstable and collapse. 
Our results suggest instead that turbulence alone is insuf- 
flcient for producing a hiera rchy of s tructu res whose mean 
densities scale with size in a iLarsonl ([198l|)-like way (recall 
our simulations are constructed to have mean densities in- 
versely proportional to their sizes, but turbulence alone ap- 
pears incapable of producing regions like run Ms8J2 within 
run Ms24J6). Turbulence must be considered in conjunction 
with self-gravity for the production of the denser structures. 

This fact is also likely to have an implication for the 
interp retation of the recent findings by Ijoung fc Mac Low! 
(J2006t ). These authors concluded from non-self-gravitating 
simulations of the supernova-driven ISM with a fixed im- 
posed supernova rate, that this driving alone is insufficient 
to sustain itself, as it does not deposit a sufficiently large 
amount of mass in Jeans-unstable regions to be consistent 
with the imposed supernova rate. Our result that gravity 
participates in the formation of Jeans-unstable regions sug- 
gests that the discrepancy between the applied supernova 
rate and the rate of production of Jeans-unstable regions in 
their simulations they reported (roughly an order of mag- 
nitude) may actually be an upper limit, so that supernova 
driving may be more efficient in driving secondary star for- 
mation than they concluded. 

Finally, the fact that the density PDF develops a high- 
density tail in the presence of self-gravity also has impli- 
cations. Several models and calculations (e.g., Elmcareenl 
I2OO2I ; [Padoan fc Nordlund 2002; Krumholz fc McKee 200^ 
rely on this distribution, and it is necessary to assess the 
degree to which such a deviation may alter the results of 
these models. 



6.4 Star formation rate per free-fall time 

We have also measured the SFRff in our three simulations, 
all of which have turbulent driving applied at the largest 
scales in each simulation, and approximately the same virial 
parameter a, thus matching the assumption made by KM05 
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to derive their equation (30) (eq. [5] in the present paper) 
that star-forming clouds and clumps are all nearly virial- 
ized, with a ~ 1 . We found the SFRfis in our simula- 
tions to be marginally consistent (i.e., within the range 
allowed by our uncertainties) with the prediction of the 
model. This is not surprising, since our suite of simula- 
tions was specifically designed to satisfy the q ~ 1 as- 
sumption of KM05. Note however that, in spite of this spe- 
cific design of the simulations, at face value the slope of 
the SFRfi vs. Ms we find is larger than the prediction by 
KM05. This may be a reflection of our result that the 
sub-boxes of the large-scale simulation Ms24J6 with the 
same size and mean density as the small-scale simulation 
Ms8J2 difi'ered from it in that they contain a non-zero mean 
convergence of the velocity field. This result contradicts 
the hypothesis of the KM05 model (and man y others, e.g., 
Elmcgrccn 2002; McKcc & Tan 2003; Mac L ow fc KlessenI 
2004 iHennebelle fc Chabrieij bOOS ) that the turbulence at 
every scale is completely random, and therefore acts as an 
isotropic pressure that provides support against the self- 
gravity of the entire structure. If in reality clouds and clumps 
contain a significant amout of kinetic energy in convergent 
motions that do not oppose gravity, the dependence of the 
SFE on Ma ai a ~ 1 might actually be somewhat steeper 
than the prediction by KM05. This in fact may be the reason 
why our simulations give a mean larger exponent in the rela- 
tion between SFRfj and Ma, since, even though each one of 
our simulations was constructed as to have a random global 
turbulent velocity field, the clumps formed self-consistently 
within them do not, on average. 

A final note of caution, however, is that our results 
have been obtained in the simplest possible numerical 
setup, namely non-magnetic, isothermal media without stel- 
lar feedback, and as such, cannot be considered definitive. 
We plan to investigate thermally bistable, magnetized media 
with stellar feeedback in future works. 
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